Nonconcave entropies in multifractals and the thermodynamic formalism 
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We discuss a subtlety involved in the calculation of multifractal spectra when these are expressed 
as Legendre-Fenchel transforms of functions analogous to free energy functions. We show that the 
Legendre-Fenchel transform of a free energy function yields the correct multifractal spectrum only 
when the latter is wholly concave. If the spectrum has no definite concavity, then the transform 
yields the concave envelope of the spectrum rather than the spectrum itself. Some mathematical 
and physical examples are given to illustrate this result, which lies at the root of the nonequivalence 
of the microcanonical and canonical ensembles. On a more positive note, we also show that the 
impossibility of expressing nonconcave multifractal spectra through Legendre-Fenchel transforms of 
free energies can be circumvented with the help of a generalized free energy function, which relates 
to a recently introduced generalized canonical ensemble. Analogies with the calculation of rate 
functions in large deviation theory are finally discussed. 

PACS numbers: 05.45.Df, 64.60.Ak, 65.40.Gr 
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I. INTRODUCTION 

Invariant measures generated by nonlinear and com- 
plex dynamical systems often show striking scaling and 
self-similar features that are reminiscent of fractals. How- 
ever, contrary to ordinary fractals, whose geometric 
structure is characterized by a single number (the fractal 
or Hausdorff dimension the scaling and self-similar 
properties of measures are usually not captured by a sin- 
gle dimension, say a, but by an infinite set of fractal or 
singularity dimensions that defines the so-called spectrum 
of singularities /(a), also known as the multifractal spec- 
trum [2, 3]. The word "multifractal" has been coined [H 
in this context precisely to suggest that a measure hav- 
ing multiscaling properties can be pictured abstractly as 
a superposition of many "pure" fractals, each having a 
dimension a and a corresponding "weight" f(a) in the 
superposition. 

To be more specific, consider a measure fi defined on a 
rf-dimensional space X. Generalizing the approach fol- 
lowed in fractal geometry, we proceed to partition or 
"coarse-grain" the space X in small boxes of equal size e 
and volume e d . The measure contained in each box is 



Pe 



(1) 



box 



and from this quantity, a local fractal dimension on, also 
called a crowding index, is defined by using the fact that 
p e% i is expected to scale as p e ^ ~ e ai in the limit where 
the boxes' size e goes to zero. Now, to account for the 
fact that cti is not constant over the partition but varies 
in general from one box to another, we count the number 
n e (a) of boxes in the partition whose local dimension is 
equal to a. From n E (a), the multifractal spectrum / (a) is 
then simply defined through another scaling relationship, 
namely n e (a) ~ e~f( a > as e — > 0. 

The multifractal spectrum f(a) is not a quantity which 
is easily calculated analytically or numerically, since it 
requires the enumeration of all the boxes in the partition 



of X having a crowding index a lying in some interval 
[a, a + Aa\. A more manageable quantity which can be 
related to f(a) is the so-called free energy function r(q) 
defined by the scaling relationship Z s (q) ~ e T ( q \ e — ► 0, 
where 



<!<\ j 



(2) 



is the partition function associated with the partition X 
of (i (the sum above runs over all the boxes of the par- 
tition with p e j ^ since q can be negative). The calcu- 
lation of r(q) parallels the calculation of free energies in 
statistical mechanics in that, if f(a) is known, then r(q) 
can be calculated as the Legendre-Fenchel (LF) trans- 
form of f(a) 0; in symbols, 



t(q) = inf {qa 



/(«)}. 



(3) 



The result that we shall study in this paper is the in- 
verse result, namely that if r(q) is known, then /(a) can 
be calculated from r(q) by taking the LF transform of 
the latter function; in symbols, 



/(«) 



inf {qa 



r{q)}- 



(4) 



This result first appeared in Refs. 0, IE], and has been 
used extensively since then to calculate the multifrac- 
tal spectrum of many phenomena, including turbulence 
H, U H, H, OH HH , geophysical processes, such as cloud 
formation and rain precipitations JSii, and fluc- 
tuations in financial time series [3j, [15( , among many oth- 
ers [l6|]. Unfortunately, there is one aspect of Eq.([3]) 
which is often overlooked when deriving it and applying 
it, namely that it can only produce concave multifractal 
spectra, since LF transforms can only yield concave func- 
tions. This basic property of LF transforms does not 
affect, as such, the calculation of r(q) from f(a) because 
it can be proved that r(q) is an always concave function 
of q. For calculating the multifractal spectrum, however, 
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FIG. 1: (a) Multifractal spectrum f(a) for the invariant density of the Ulam map. (b) Corresponding free energy function 
r(q). (c) Legendre-Fenchel transform of r(q). 



there is a problem because f(a) need not be concave, 
which means that /(a) cannot always be calculated as 
the LF transform of r(q). 

Our goal here is to illustrate these observations with a 
number of basic examples, and to state the precise con- 
ditions, based on convex analysis, that ensure that f{a) 
can be calculated as the LF transform of r(q). These con- 
ditions will be discussed in the context of four physically- 
relevant multifractal models: one related to turbulence, 
another related to diffusion- limited aggregates, and two 
others related to chaotic systems. In an attempt to offer 
a workable solution to the problem of calculating noncon- 
cave multifractal spectra, we shall also study a recently 
introduced generalized canonical ensemble to show that 
nonconcave spectra can be obtained from a modified ver- 
sion of the LF transform. This part will actually provide 
an explicit calculation of a nonconcave spectrum f(a) 
which uses this modified LF transform. We shall com- 
ment finally, in the concluding section of the paper, on 
analogies between nonconcave multifractal spectra, non- 
concave entropies in statistical mechanics, and large de- 
viation theory. 



II. TWO SIMPLE EXAMPLES 

We begin by considering two explicit examples of mea- 
sures whose multifractal spectra are not given by LF 
transforms of their free energy functions. The first exam- 
ple was previously discussed in Ref. [ljj (see also Ref. 0), 
and will serve here as a starting point to our discussion of 
the validity of the LF transform of ^j. The measure or, 
rather, the density in this case that we consider is given 
by 



p(x) 



71- VT 



~2' 



(5) 



where x S [—1,1]. This density arises as the invariant 
density of the Ulam map and the Tchebyscheff maps. 
Applying a partition of size e on the interval [—1,1], it 
can be seen that the two boxes of the partition located 
near the boundary points x — ±1 have measure p e% i — 
ep(x) — e 1/2 , so that f(a) = at a = 1/2. All the other 
boxes have measure p £ ^ ~ e, so that /(a) = 1 at a = 1, 



as there are approximately n e (a) ~ e~ x of these boxes. 
Combining the two results, and setting n e (a) = for 
a $ {1/2, 1}, we obtain 



/(«) = 



-co 



a = 1/2 
a = 1 
otherwise. 



(6) 



This spectrum is shown in Fig. [TJa). 

At this point, we go on to prove that f(a) cannot be 
expressed as the LF transform of r(q) by direct calcula- 
tion. Starting from the asymptotic (e — * 0) expression of 
the partition function 



Z E (q)=Y,Pli~ £q/2 + £ ~ l£q > 



(7) 



we first find 



T(q) = mm{q-l, q /2} = { q q /2 1 «>2 (g) 

Then keeping track of the two separate regions q > 2 and 
q < 2, we find 

inf {qa - r(q)} = inf ( q< f ~ f ) q J I 

<jsr q£R { q(a - 1) + 1 q < 2 

2a -I ae [1/2, 1] 
— oo otherwise. 

Comparing this result with Eq.©, we see that f(a) cor- 
responds to the LF transform of r(q) for a $ (1/2,1) 
only; see Figs. [T^a) and^c). For a € (1/2,1), the LF 
transform of r(q) is finite, while the true spectrum /(a) is 
formally equal to — oo, as there is no box in the partition 
of p(x) with local exponent in the range (1/2, 1). 

This example can be generalized to illustrate another 
problem when trying to obtain f(a) from r(q). Consider 
a dynamical system in (i-dimensions whose invariant den- 
s' 4 , except at a 



sity is everywhere finite, so that pi jE 
of singular points whe 
£k < L The partition function for this 



finite number k of singular points where pi :E ~ e Qi = e ^ 



with £i,£ 2 , 
density is 

Ze(q) = ] 



.Pi 



-d^dq_ 



-e 



(10) 
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so that 



r(q) = min{(q - l)d, d^q, d£ 2 <7, d£kO}- 



(11) 



The minimum can be calculated explicitly and yields 



(d(q-l) g<(l-r) _1 
TW -\d?q q>(l-C)-\ 



(12) 



where £* = min^. We see here that the function r(q) 
"overlooks" all the singularities £j , except for the smallest 
one. Therefore, any perturbation of the singularities £j 
that keeps £* invariant will change f(a) but not r(q), 
which implies that the mapping of r(q) to f(a) must be 
indeterminate as there is an infinite number of spectra 
associated with the same free energy. Physically, this 
also implies that r(q) does not offer the most complete 
description of the dynamical system, since this function 
overlooks, as we said, all but one singularity. To really 
obtain a complete picture of all the singularities of the 
system, one must resort to calculate f(a) and not just 
r(q)- 



III. THEORY OF LF TRANSFORMS 

The results of the two previous examples are very sim- 
ple and show at once that f{a) cannot in general be 
expressed as the LF transform of r(q), contrary to what 
is claimed in most if not all references on the subject. 
The problem, as was mentioned, is that LF transforms 
can only yield concave functions, which means that these 
transforms cannot be used to calculate nonconcave multi- 
fractal spectra, including those of the two examples con- 
sidered before. To make this observation more rigorous, 
we introduce in this section a few concepts and results of 
convex analysis, beginning with the concept of support- 
ing lines. (All the definitions and theorems discussed here 
can be found in Ref . [l8[ ; see also Chapter VI of Ref . (l9j 
and Appendix A of Ref. (2fjj.) 

Definition 1. A function / : R — > K admits a supporting 
line at a if there exists a constant n such that 



f(P)<f{<x)+ri[p-a) 



(13) 



for all (3 G 



This definition means graphically that we can draw a 
line on top of the graph of f(a) that does not go under 
that graph (see Fig. (2); hence the word "supporting." 
With this picture in mind, it is easily seen that, if / 
admits a supporting line at a and is differentiable at a, 
then the slope 77 of the supporting line must be such that 

/'(«) = * 

The importance of supporting lines comes from their 
association with LF transforms, and from the fact, more 
precisely, that they determine whether such transforms 
are involutive, that is, whether they are their own inverse. 
In the context of f{a) and T(q), this means precisely the 
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FIG. 2: (a) A generic nonconcave multifractal spectrum f(a) 
(full line) together with its concave envelope f**(a) (dashed 
line). The two functions coincide outside the open interval 
(a,i,ah)- The point a of the multifractal spectrum admits a 
supporting line (concave point), while the point b does not 
(nonconcave point). 



following. First, recall that r(q) can always be expressed 
as the LF transform of f(a), so Eq.([3]) is always valid in- 
dependently of the shape of f(a). This follows essentially 
from the fact that r(q) is an always concave function of 
q @, a fact that can be proved using Holder's inequality. 
The inverse transform shown in 0, however, is not gen- 
erally valid, and this is where supporting lines become 
important, as expressed in the next theorem. 

Theorem 2. If / admits a supporting line at a, then 
/ at a can expressed as the LF transform of r(g) as in 
Eq.Q. In this case, we say that / is concave at a. On 
the other hand, if / does not admit a supporting line at 
a, then / at a does not equal the LF transform of r(q). 
In this case, we say that / is nonconcave at a. 

The two complementary results expressed in the the- 
orem above are usually rephrased in convex analysis by 
defining the function 



f**(a) = mi{qa - r(q)}. 



(14) 



In terms of f**(a), we then have the following result. 

Theorem 3. f(a) = f**(a) if and only if / admits a 
supporting line at a. 

For the remaining, it is useful to note that /**(«) cor- 
responds in general to the smallest concave function sat- 
isfying f(a) < f**(a) for all a e M. For this reason, 
f**(a) is called the concave envelope or concave hull of 
f(a). This implies, in particular, that if f(a) admits no 
supporting lines over some open interval, say (07,0^) as 
in Fig. [21 then f**(a) must be affine over that interval, 
by which we mean that f**(a) has a constant slope over 
that interval. This last property, which is related to the 
Maxwell construction [2l|, [22|, is illustrated in Fig. [2j 

All of the properties of /(a) and f**(a) in relation to 
LF transforms can be verified for the two examples con- 
sidered previously. In the case of the invariant density 
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of the Ulam map, for example, the concave hull of f(a) 
is the function displayed in ©; it is obviously such that 
f{o) < f**( a ) an d is concave contrary to f(a). More- 
over, it is easily verified from Fig. [1] that the two points 
a = 1/2 and a = 1 admit a supporting line, which ex- 
plains why /(a) = f**(a) there. These two points admit 
in fact an infinite number of supporting lines. For the 
point a — 1/2, for example, all lines attached to (1/2,0) 
with slope in the interval [2, oo) are supporting in the 
sense of (fl~3|) . For a = 1, the supporting lines have slopes 
in the interval (— oo, 2]. 

We can go further in our analysis of f(a) and r(q) by 
calling attention to the fact that 

r(q)= M{qa-f**(a)}. (15) 

Therefore, r(q) is not only the LF transform of f(a), as 
stated in Eq.([3]), but also the LF transform of f**(a). 
This result is general: it holds for any function f{a) and 
its concave envelope f**(a) defined as in Eq. (fl4|) as the 
double LF transform of f{a) or, more compactly, as 

f*=r* = (/*)*, (16) 

where the star stands for the LF transformation. To 
summarize, we then have r = /* , r* = (/*)* = /** and 
(t*)* = (/**)* = /* = r. This chain of equalities can be 
expressed in a more transparent way using the following 
diagram: 

f(a) Ar(^f(a), (17) 

which makes obvious the fact that there may be more 
than one spectrum related the same free energy. In fact, 
all f(a) having the same concave envelope lead to the 
same r(q), as can be verified in the second example con- 
sidered before. Finally, note that the chain of equalities 
reduces to r = /* and t* = /, or equivalently to 

r(q) ± f(a), (18) 

when /(a) = /**(a) for all a e K, that is, when /(a) is 
everywhere concave. 

Having listed all the relationships that exist between 
/(a), r(q) and /**(«), we can now fully address the main 
issue of this paper, which is to determine when f(a) can 
safely and completely be calculated as the LF transform 
of r(q). From the chain of equalities and diagrams shown 
above, this amounts to determine when the LF transform 
is involutive; that is to say, under which conditions does 
the diagram ((TTJ) reduce to the diagram of (fT8")l ? 

A first obvious answer to this question is given by re- 
calling what we have just mentioned about the diagram 
of l|18p. namely that if f(a) is everywhere concave, then 
the multifractal spectra f{a) can completely be calculated 
as the LF transform of the free energy function r(q). As 
such, this answer is complete but not very practical be- 
cause it is based on f(a) and so presupposes that we 
know f(a). A more useful criterion can be stated from 
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FIG. 3: Free energy function r(q) associated with the mul- 
tifractal spectrum f(a) shown in Fig. [2] The LF transform 
of the concave envelope f**(a) of /(a) yields the same free 
energy function. 



the point of view of r(q) alone by using a result of con- 
vex analysis connecting nonconcave or affine regions of 
f(a) with nondifferentiable points of r(q). This result is 
stated next without a proof; see [IH for more details (see 
especially Theorems 23.5 and 26.3). 

Theorem 4. Suppose that f(a) is nonconcave over some 
open interval (ai,ah) (Fig. ^ or that f(a) is concave 
but affine over (pu, ah). Then r(q) is nondifferentiable at 
some critical value q c corresponding to the slope of /** (a) 
over the interval (a/,a^). Moreover, the left- and right- 
derivatives of r(q) at q c equal ah and a;, respectively 
(Fig.©. 

From this result, we arrive at our criterion by taking 
the contrapositive: if r(q) is everywhere differentiable, 
then f(a) is concave everywhere with no affine parts. 
Thus, from the point of view of r(q), f{a)can completely 
be calculated as the LF transform of r(q) if the latter 
function is everywhere differentiable. Taking the view 
that nondifferentiable points of r(g) represent first-order 
phase transitions for multifractals [23|, [24], [2f| [111 l27l ]. 
this is equivalent to saying that f(a) can completely be 
calculated as the LF transform of r(g) in the absence of 
first- order phase transitions. If there is a first-order phase 
transition, then either f(a) is nonconcave somewhere, in 
which case / r* , or else f(a) is affine somewhere, 
in which case / = r*. Unfortunately — and this is an 
important point — there is no way to distinguish the two 
cases from the sole knowledge of r(q) (see Figs. [2] and [3]). 
Thus, if r(q) has one or more nondifferentiable points 
and if there is no reason to think that f(a) is concave, 
then one must resort to calculate f(a) by means which 
do not rely on r(q). 



IV. APPLICATIONS 

We now revisit some examples of multifractal models 
that have been discussed in the physics literature, and 
point out where and how our results of the previous sec- 
tion apply. 



5 



A. Multifractal turbulence models 

One of the first field of study for which multifractal 
models have been developed is fully developed hydro- 
dynamic turbulence [J, 0]. The basis of these mod- 
els is that, in the turbulent flow, velocity increments 
5v(l) = \v(x + I) — v(x)\ at a given distance / scale with 
local exponents h, which are distributed on a fractal set 
with fractal dimension D{h). This notation is taken from 
Benzi et al. 0] , and can be translated to our notation us- 
ing the following identifications: 

h = a 

I = e 

3-D(h) = -/(a) 

p = q 

Cp = r(q). (19) 

Here ( p denote the scaling exponents of moments of ve- 
locity increments in the inertial range, 

((5v) p ) ~ fi*. (20) 

In multifractal turbulence models, the probability to ob- 
serve a local exponent h is given by 

Pl(h) ~ l 3 - D (' l \ (21) 

which is equivalent to our notation n E (a) ~ e^^ a \ 
Moreover, for the scaling exponents p one has 

C P = min{/ip + 3 - D(h)}, (22) 

h 

which is equivalent to the LF transform of Eq.Q. In 
practice, one extracts D(h) from the scaling exponents 
£ p , which can be measured in experiments. In our no- 
tation, this means that one determines f(a) by the LF 
transform of the experimentally measured r(g). 

In view of all the results derived before, we can notice 
that multifractal turbulence models in their current form 
can only deal with the convex hull of the spectrum of 
singularities. The true spectrum D(h) of a turbulent flow 
is fully determined by the underlying dynamics, i.e., the 
Navier-Stokes equation, and there is a priori no reason 
to think that this spectrum should be a concave function 
of h. Our arguments of the previous section now allow 
for an experimental check of the concavity of D(h): if 
the experimentally-measured £ p is differentiable within 
the precision allowed by the experiment, then D(h) is 
concave and therefore given by the LF transform of (, p . If 
one observes that Cp is nondifferentiable, then either D(h) 
is nonconcave or else is affine. Both cases are consistent 
with the fact that £ p is nondifferentiable, but there is no 
way to tell from £ p which one of the two spectra is the 
actual one. 

The exponents £ p have been measured in many exper- 
iments; see, e.g., Refs. d, H[. Within the experimental 
uncertainties, they are usually described by a smooth 



function of p, although one cannot fully exclude the exis- 
tence of phase transitions. The data found in Refs. [!,[§], 
for example, show a relatively strong change of slope near 
p = 3. If the existence of such transitions were confirmed 
(e.g., via the study of theoretical models of turbulence), 
then one would have to check that the underlying spec- 
trum D(h) is either nonconcave or concave but affine 
somewhere. If D(h) is affine (see, e.g., Ref. 28]), then 
that spectrum is correctly given by the LF transform of 
Q p . If D(h) is nonconcave, then that spectrum is not 
fully given by the LF transform of C, p . In this case, ( p 
cannot provide a complete description of the turbulent 
flow, since all the local exponents h in the nonconcave 
region of D(h) are "overlooked," in the spirit of the two 
simple examples discussed before, by Cp- I n short, these 
exponents are described by D(h) but not by £ p . 



B. Diffusion-limited aggregates 

For the second example, we consider multifractals as 
generated by diffusion-limited aggregates (DLA) [2^, [3(| 
[3lj ]. Jensen et al. [3l[ provide convincing evidence that 
the function r(q) calculated for the harmonic measure of 
their DLA cluster exhibits a first-order phase transition 
at q = — 0.23± 0.05 (see their Fig. 3). From this result, 
all that can be said about their /(a) spectrum obtained 
by LF-transforming r(q) (shown in their Fig. 4) is that 
it is the concave hull f**(a) of the true f(a) spectrum. 
Note that the spectrum displayed in Fig. 4 of Ref. [3l| 
does show an affine part, so it is consistent with the fact 
that r(q) has a nondifferentiable point. However, there 
is a priori no reason why DLA clusters should possess a 
concave spectrum of singularities, so that the part where 
f(a) is seen to be affine could just as well be nonconcave. 
Therefore, at this point we may conclude that the true 
/(a) spectrum of the DLA cluster studied by Jensen et 
al. is as yet unknown. 

It should be remarked that, although Jensen et al. pro- 
vide evidence to the effect that r(q) possesses a nondif- 
ferentiable point, the r(q) which is calculated in practice 
is actually always analytic if one deals with finite-size 
DLA clusters. The spectrum f**(a) which is calculated 
from r(q) is, in this case, necessarily concave and has no 
affine parts. The nondifferentiable point of r(q) and the 
concomitant affine part of f**(a) appear, formally speak- 
ing, only in the "thermodynamic" limit of infinitely large 
clusters. 



C. Chaotic systems 

Multifractal spectra have been calculated for many 
examples of chaotic dynamical s yst ems, including the 
Henon map [13] (see also Ref. [33j]), and the driven 
damped pendulum [34]. The particularity of these two 
examples is that they seem to give rise to first-order phase 
transitions, which are also referred to as (/-phase tran- 
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sitions. Accordingly, the question arises as to whether 
these phase transitions emerge out of a nonconcave f(a) 
or an affine /(a). 

The question, as it stands, is not resolved in the pa- 
pers that treat these examples because they assume that 
f(a) is always the LF transform of r(q), which means 
that they implicitly assume that f{a) is always concave. 
It must be observed that some of the reported spectra ap- 
pear to be affine (see, e.g., Fig. 2 in Ref. [32[), so they are 
not problematic — they satisfy the concavity assumption. 
However, some other spectra are clearly nonconcave; see, 
e.g., Fig. 1 in Ref. [HI and Figs. 2 and 3 in Ref. [H. For 
these, we must be careful because the observed noncon- 
cavity could be a finite-size effect inherent to the fact that 
f(a) is computed numerically for a finite coarse-graining 
resolution s. Thus it could be that the nonconcavity of 
f(a) observed for e > disappears as e — ■> 0. To verify 
this, one would need to perform a finite-size analysis of 
the data by computing /(a) for decreasing values of e and 
study the convergence of the results. Similar finite-size 
analyses have been performed in the context of numerical 
calculations of the microcanonical entropy function; see 

Refs. mm 13. 



D. Spectrum of dynamical indices 

First-order phase transitions have also been studied in 
the context of chaotic systems at the level of an entropy- 
like quantity referred to as the spectrum of dynamical 
indices or expansion-rate spectrum 0, [H, US HI] • These 
phase transitions typically occur for nonhyperbolic dy- 
namical systems. 

It would take us too far to explain the notion of dy- 
namical indices and the many examples for which this 
quantity has been studied. Let us only mention that the 
spectrum of dynamical indices is a dynamical analog of 
/(a), and that affine and nonconcave spectra of dynami- 
cal indices have been reported in the literature; see, e.g., 
Refs. m Him EUS El- Most of these references, un- 
fortunately, share the same problem as those discussed so 
far: they assume that the spectrum of dynamical indices 
can always be calculated as the LF transform of a dynam- 
ical analogue of the free energy function r(q). Most of 
them assume this even when reporting the computation 
of nonconcave spectra; see Yoshida and Miyazaki [43[ for 
a noticeable exception. 



generalized form of free energy function. The method 
was proposed recently in Refs. [2(| EH in the context of 
nonconcave entropies within the microcanonical ensem- 
ble, and will be illustrated here in the context of the first 
example considered in Sec. [ill 

The form of generalized free energy that we shall con- 
sider is based on a generalization of the partition function 
given by 



Ze{q,g) 



£ qai+g{oti) 



(23) 



where on represents the local fractal exponent associated 
with the probability p e ,i, and g is an arbitrary smooth 
function. This new form generalizes the standard canon- 
ical partition function, in the sense that 



(24) 



For definiteness, we shall adopt the choice g(a) = 7a 2 
with 7 € K. Therefore, the generalized partition function 
that wc consider is 



(25) 



We call this partition function the Gaussian partition 
function; its corresponding free energy function 



lim 

e— 



lnZ e (g,7) 
lne 



(26) 



is called the Gaussian free energy. Note that this new 
free energy is a function of two real parameters, q and 7, 
and that r(q, 7 = 0) = r(q). 

The rationale for generalizing the standard free energy 
function r(q) to r(g, 7) is that it modifies the structure 
of the LF transform which connects r(q) with /(a), and 
thus modifies the conditions which ensure that f(a) can 
be written as the LF transform of a free energy function. 
We spare the reader with the details of this modification 
which can be found in Refs. [2(||44j|. F° r our purpose, we 
shall only note the generalized versions of the LF trans- 
forms that connect r(q, 7) and f(a); they are given by 



r(q, 7) = inf{qo! + 7a 2 - f(a)} 



(27) 



and 



V. GENERALIZED FREE ENERGY 
FUNCTIONS 

At this point, we have emphasized more than once that 
a nonconcave multifractal spectrum cannot be calculated 
as the LF transform of its corresponding free energy func- 
tion. Our goal in this section is to offer a practical so- 
lution to this problem by illustrating a method for ob- 
taining nonconcave spectra through LF transforms of a 



f(a) = M{qa 

1:1 



- 7« 



"(9,7)}. 



(28) 



The first LF transform holds, like its standard version 
(7 = 0), for any spectrum /(a), be it concave or not. 
The surprising virtue of the second LF transform is that 
it also holds true for basically any /(a), contrary to the 
standard version (7 = 0) which applies only when f(ct) is 
concave. Rather than proving this result, we shall verify 
that it is valid for the nonconcave multifractal spectrum 
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shown in Fig. QJa). That is, we shall obtain that noncon- 
cave spectrum by inverting, in the manner of Eq. (|28p , its 
associated Gaussian free energy r(q,j). 

First, we calculate r(q, 7) starting from Z e (q r , 7): 



«/2+7/4 



e -i e? +7_ (29) 



Taking the limit e — > yields 

r(g, 7) =min{ (Z /2 + 7/4,9 + 7-1}. (30) 

The solution of the minimum can be found explicitly; it 
has the form 



■^-{llV-'t 



q > <? 7 

q < g 7 , 



(31) 



where q 7 = —37/2 + 2. Next, we apply formula 
using this solution for r(q, 7 ). This leads us to solving 
the following variational problem: 



qa + 7a 2 — q/2 - 7/4 q > g 7 



7 = inf \ 



7+1 g < ? 7 
Grouping the variables together, this is equivalent to 
g(a-l/2)+ 7 (a 2 -l/4) q > q y 



(32) 



g(a-l) + 7 (a 2 -l) + l 



q<q y 



(33) 



Let f\(a,q, 7) denote the top expression in the brackets 
and f2(a,q,"f) the lower one. With this notation, it can 
be noted that, for a = 1/2, 

/ 2 (l/2,g, 7 ) >/i(l/2, 9)7 ) = (34) 

for all q < q 1 and 7 € R. Therefore, 

inf / 2 (a,g,7) = /i = ) (35) 

9<9 T ,7 

and 7 = at a = 1/2. Similarly, for a = 1, we have 

/i(l,g J 7)>/a(l,g,7) = l (36) 

for all <z > 97 and 7 6 I, so that 7 = / 2 = 1 at a = 1. 
For all other values of a, it is possible to set q and 7 in 
such a way that I = —00. At the end, we are left with 



a = 1/2 
/ = <( 1 a = 1 

-00 otherwise, 



(37) 



which is the precise expression of /(a), as given in Eq.©. 
Therefore, we have shown that this nonconcave spectrum 
can be expressed as in Eq. I|28p as a modified LF transform 
of a generalized free energy. 

The same method can be applied to calculate other 
nonconcave spectra. In fact, it has been conjectured 
that the method can be used to calculate any noncon- 
cave spectrum (viz., nonconcave entropy function) as the 
LF transform of a properly-chosen generalized free en- 
ergy. More details about this universality property can 
be found in Refs. [II EH- 



VI. CONCLUSION 

We have shown in this paper that one must be careful 
when calculating the singularity spectrum of multifrac- 
tals as the LF transform of its corresponding free en- 
ergy, since LF transforms can only yield concave func- 
tions. This word of caution has implications for most of 
the studies published so far on multifractals, including 
those on multifractal models of turbulence, as they have 
taken for granted that the multifractal spectrum is the 
LF transform of the free energy no matter what the spec- 
trum looks like. This, as we have seen, is only true if the 
spectrum is concave; if it is nonconcave, then one must 
resort to calculate it directly from its definition. Another 
possibility is to use a generalization of the canonical en- 
semble which can be used to extract nonconcave entropies 
from a generalized version of the free energy function. 
This way of doing was sketched here in the context of 
a simple example of nonconcave multifractal spectrum, 
and is presented in full details in Refs. [20L PH] . 

In the end, it should be noted that the results that 
we have discussed in this paper are not special to mul- 
tifractals, but apply actually to any field of investiga- 
tion which uses LF transforms. In statistical mechanics, 
for example, the LF transform that connects the entropy 
function of the microcanonical ensemble with the free 
energy function of the canonical ensemble becomes non- 
involutive when the entropy is nonconcave. When this 
happens, we say that there is nonequivalence of ensem- 
bles [i^, [4(| , since one is then unable to obtain the true 
entropy function of the microcanonical ensemble solely 
from the knowledge of the free energy of the canonical 
ensemble. The notion of generalized canonical ensemble 
has been developed precisely in this context. 

Similarly, in large deviation theory, it has been known 
for some time that nonconvex rate functions cannot be 
calculated by means of LF transforms of functions analo- 
gous to free energy functions [13, IH, [4t| [5(| • A multifrac- 
tal spectrum is in essence an entropy function, and an en- 
tropy function is in essence a rate function [5ll[52l [53l.[54|. 
so there is actually a deep connection with what we have 
presented here and what is known in large deviation the- 
ory [50l | . For example, the result relating the differen- 
tiability of r(q) and the (strict) concavity of f(a) can 
be put in correspondence with a result of lar ge d eviation 
theory known as the Gartner-Ellis Theorem [l9l l49l [55| . 
Furthermore, the result stating that the LF transform 
of f(a) always yields r(q) (see Sec. IIII|) can be put in 
correspondence with a large deviation result known as 
Varadhan Theorem; see Refs. (l9ll55|. 

From these correspondences, it is but a small step 
to conjecture that nonconcave entropies or, more gen- 
erally, nonconvex rate functions should show up in other 
physical theories in which large deviations are at play. 
One such theory that comes to mind is the thermody- 
namic formalism of dynamical systems 0, [56j]; another 
is nonequilibrium statistical mechanics. We have already 
alluded to the first theory when discussing the spectrum 
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of dynamical indices in Sec. IIVI Concerning nonequi- 
librium statistical mechanics, the reader will find an ex- 
ample of nonconvex rate function in a recent paper by 
Imparato and Peliti [57| (see their Fig. 12). This paper 
discusses thermodynamic fluctuations in systems driven 
out of equilibrium. A large deviation result in that con- 
text is referred to as a fluctuation theorem, whereas a 
rate function is called a fluctuation function. 



one by Riedi [60( in which examples of nonconcave mul- 
tifractal spectra are discussed. 
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